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Kinetic pinning and biological antifreezes 
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Biological antifreezes protect cold-water organisms from freezing. An example are the antifreeze 
proteins (AFP's) that attach to the surface of ice crystals and arrest growth. The mechanism for 
growth arrest has not been heretofore understood in a quantitative way. We present a complete 
theory based on a kinetic model. We use the 'stones on a pillow' picture Our theory of the 

suppression of the freezing point as a function of the concentration of the AFP is quantitatively 
accurate. It gives a correct description of the dependence of the freezing point suppression on the 
geometry of the protein, and might lead to advances in design of synthetic AFPs. 

PACS numbers: 87.14.Ee, 68.08.-p, 81.10.Aj 



In polar regions, many fish, insects and plants flour- 
ish at temperatures well below the freezing point of their 
bodily fluids 0, Q. Often, particularly in fish this 
is a non-colligative effect, namely that ice crystal growth 
within the organisms is arrested by a class of plasma pro- 
teins called anti- freeze proteins (AFP). They are usually 
peptides or glycopeptides. These molecules attach irre- 
versibly to ice surfaces |5j and lead to arrest of crystal 
growth^, Q until the water containing the AFP's is su- 
percooled by as much as 2 degrees C. The mechanism for 
the suppression of the freezing point is not understood. 
In this paper we present a growth model allows us to un- 
derstand AFP's in considerable detail, and, in particular 
to calculate the dependence of the undercooling on AFP 
concentration. 

Growth arrest by AFP's occurs because the protein 
adsorbs on the surface of the growing ice crystal, is not 
incorporated in it and suppresses growth at that site. 
One version of this notion has been called the 'stones 
on a pillow model' 0, 0] • It assumes a thermally rough 
crystal, appropriate for most of the surface of ice near 
its freezing point. (The basal plane of ice is facetted, 
but AFP's adsorb mostly on the prism planes @.) The 
AFP's are obstacles to growth so that the crystal surface 
bulges between the attachment sites, and the freezing 
point is depressed by the well-known Gibbs-Thompson 
effect, namely that a curved surface has a lower freezing 
point than a flat one; see Figure (la). This is the physical 
picture we will adopt. The surface of the crystal which is 
not under the AFP must have constant mean curvature, 
k. The Gibbs-Thompson condition [2j is: k = —5t/Io 
where St = (T m — T)/T m (the undercooling) and T m is 
the equilibrium melting temperature [i|: Also, l = 7/ A 
is a characteristic length, 7 the interfacial tension, and 
A is the latent heat of fusion per unit volume. For ice 
l Q « 1A. 

Consider a surface, h(x, y) which has stopped grow- 
ing. We model the AFP's as spheres of radius b ~ 10A. 
(We treat more complicated molecular geometries, be- 



low.) Assume that the particles are half buried in the 
ice, and fix the ice surface at the equatorial plane of the 
sphere. The angle of the ice with the equatorial plane can 
take on any value. We use the small slope approximation, 
\dh/dx\ « 1, \dh/dy\ << 1, so that the curvature of the 
interface is given byK~ Ah. Then the Gibbs-Thompson 
condition becomes the Poisson equation: 



Ah = -5 T /l = -An p. 



(1) 



This equation is familiar in electrostatics, h plays the 
role of the potential, and the effective (positive) charge 
density, p is defined by Eq. (JJJ. 

For simplicity, consider a periodic square array of 
AFP's. The boundary condition on h at the edge of the 
unit cell is that the normal derivative vanishes, and at the 
edge of the AFP, h = 0. This problem is easily solved 
numerically and a contour plot is shown in Figure (lb). 
In order for a solution to the equation to exist, the AFP 
must act as a negative charge so that the system is neu- 
tral. That is, p = —nq, where n is the density of AFP's 
on the surface. By integrating Eq. over the surface 
we find: 



((r -!•(). V/i) ]p _ r4 | =8 T /(2nl n) 



(2) 



where (•) is the average around the edge of the AFP. The 
effective charge on the AFP is: 



g = -(l/2)((r-r i )-Vft), r _ r4 | 



(3) 



Note from Eq. J3J that the slope at the edge of the 
molecule increases with undercooling. It is reasonable to 
assume that if the slope exceeds some critical value, Xi 
the antifreeze molecule will be engulfed. Here x 1S se t 
by the physical chemistry of the AFP and the interface. 
We must have x ~ 1- The maximum undercooling, 5 T , is 



given by Eq. © with < dh/dr >= x- 
5 T of order 1 degree we need n = ri 
6xlO n cm~ 2 , or a distance betwee 
From Eq. © we find —2q < xb < 0. 



For b — 10A , and 
= 6* T /(2ir X l b) w 



FIG. 1: a). The 'stones on a pillow' model. AFP molecules, white circles, attach to a surface, b). Visualization of a surface of 
constant mean curvature. The size of the unit cell has been set to unity. 



Since the adsorption is thought to be irreversible, the 
relationship between n, the density of AFP on the sur- 
face, and the concentration in solution must depend on 
the kinetics of the growth of the crystal, which we now 
model. In the simplest picture, the growth speed of the 
crystal boundary is proportional to the variation of the 
overall free energy of the system with respect to the nor- 
mal displacement. In the small slope limit we take this 
to be Sh. Then: 



h = -r 



Sh 



7[l + |Vh| 2 ] 1/2 +A,5T/i 



d 2 r. 



(4) 



The first term of the integrand is the surface energy, and 
the second is the difference of the bulk free energies of 
solid and liquid phases close to melting point. T is a 
kinetic coefficient related to the rate of attachment of 
water molecules to the ice surface. 

From the calculus of variations [lfjj this equation is 
equivalent to: 



h = v a [loAh + St] 



(5) 



Here v — TA. This is a diffusion equation for h. The 
AFP's are pinning centers for the surface, whose effect 
may be expressed as a set of boundary conditions: 



\=bo 



(6) 



Here is the position of i-th AFP and hi is the height 
of its equatorial plane. 

Suppose that the AFP's in the water adsorb at random 
on the ice surface with rate, fc + , per unit area. The value 
of hi is the local position of the interface at the moment 
when the i-th AFP is adsorbed. The AFP's effectively 
disappear when they get buried under the ice surface. 

It will be useful to split h (r) into two contributions: 
h (r) = H +(f> (r), where H is the average height, and 4> (r) 
measures the small-scale variations. Suppose the AFP's 
are randomly distributed with density, n. We choose <j> 
to develop in time as: 



4> = v l (A(f> + A7rp) 



(7) 



where p — n < q > is chosen to make the surface neutral, 
and < q > is the average charge. It is not hard to show 
that the time derivative in Eq. J7J) for (j> is small (the 
quasi-static limit). Thus: 



Acj) = -Air 



P- 



(8) 



The point charges account for the boundary conditions, 
Eq. JBJ. The values of qi are not fixed: as the interface 
moves, the charges vary from at the moment of creation 
to —\b/2, at which point the AFP is engulfed by ice. 

Near an AFP, (b < |r — r;| <C lA/ra), 6 is dominated 
by the contribution of a single charge |ll| : 



<£(r) ~ -gilog7r(i 



Therefore: 



Qi 



H-hj 
log (irnb 2 ) 



Yi) n 



< 



(9) 



(10) 



Since \q\ lies between and X&/2, the 'hole' in the inter- 
face associated with the i th AFP is of order \qi\ ~ %6/4. 
Subtracting Eq, J7J) from Eq. J5J, we find: 



H = v Q [St - 47rZop] 



(11) 



Suppose the interface moves with constant speed V. 
Since each hi is fixed and H changes uniformly with 
time, Eq. IjlUfl implies that the magnitudes of individ- 
ual point charges qi are uniformly distributed between 
and — x&/2, i.e. (q) = —yhjA and p = nb/4. This implies 
that V = v (St — irxl bn). 

The evolution of n can be calculated by noting that 
its rate of increase is fc+, and that its rate of decrease 
is the rate that AFP's are engulfed. Eq. (|10|l implies 
q\ = V/ log(7rn6 2 ). Since the qi are uniformly distributed 
we must have dn/n = —qidt/(xb/2). This gives: 



2nV 



X&log (l/nnb 2 ) 



(12) 
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In the steady state we put h = in Eq. 1(12(1 . and use 
the expression for V. Thus, 

fc+ = n[2-n/n*] 
Vo^t x61og (l/7rn& 2 ) ' 

where n* = 5t (27rx&io)~ • The right hand side of this 
equation has a maximum for n ~ n*. Thus, there is no 
constant speed solution unless the right hand side is small 
enough, i.e. for small enough k + or large enough St- The 
threshold S T obeys: 

(S* T ) 2 / \og(2 X lo/bd T ) =s { X b) 2 2nl k + /v o (14) 

To test our approximations, notably the quasi-static 
limit for <j), we have performed a numerical simulation 
of Eqs. (J5J, ©, coupled to the random adsorption of 
AFP. Results on a 150 x 150 square lattice, with each cell 
representing a single AFP, are shown in Figure El These 
results support our analysis. In fact, the transition from 
the steady growth to the arrested interface regime occurs 
at a somewhat lower adsorption rate fc + than predicted. 
Since we have not performed a complete stability analysis 
of the steady-growth solution, Eq. 1(14(1 gives only an 
upper bound for the critical value of k+. However, the 
discrepancy is rather small. 

Beyond the transition point, growth stops. The result- 
ing static interface must obey Eq. © with h = 0. Once 
again, we recover the Poisson equation, Eq. (JJJ. After 
arrest, as Eq. 1)12(1 shows, n increases as irreversible ad- 
sorption continues until limited by some aspect of surface 
chemistry that we have not considered in our model. 

Real AFP's are often anisotropic. In order to ac- 
count for this, we assume that the region blocked is 
elliptical, with semi-major and semi-minor axes a and 
6, respectively. The potential h (x, y) can be found 
by conformal mapping |l2j . Using this method, it is 
easy to show that Eq. (|10|) is replaced by: H — hi ~ 

— qi/ log ( wri (a + b) 2 /4j . The slope of h will reach its 
critical value near x = ±a, and its maximum will be 
given by the radius of curvature of the ellipse, of order 
b. Thus g max = With these changes, the left-hand 

side of Eq. H14|) becomes 

(S T ) 2 /\og(8xlob/(a + b) 2 5 T ) (15) 

Eq. (|14|l involves k + /v The rate, fc+ is certainly pro- 
portional to the concentration, C, of AFP in solution: 
k + = vaC where va has the units of velocity. Therefore, 
va/v is a parameter which we need to estimate. 

In any experiment, soon after the adsorption process 
starts water near the ice surface is depleted of AFP. The 
thickness of the depletion layer depends on the diffu- 
sion coefficient D, as \/Dt. That is, fc+ becomes time- 
dependent, and thus va — y/D/t- The same argument 
applies to the rate of crystallization itself: our estimate 
for v fails as soon as the process is limited by thermal 



diffusion through a diffusive layer. Thus v a ~ yj ktT/M, 
where kt is thermal conductivity of water. From the Ein- 
stein formula D = ksT /6irr]b: 

va/v = y/lTb ~ HT 2 (16) 

where A = kB-h/6irr]KT = 1-2 X 10~ 10 cm. By combining 
Eqs. 1(14(1 and ((16(1 we can write: 

5* T [\og(2 x l /bS T )r 1/2 ~ xV^VC, (17) 

For anisotropic AFP's we must use Eq. 1(15(1 and also 
note that the longer axis, a, sets the hydrodynamic radius 
of the ellipsoidal particle. Its diffusion coefficient may be 
approximated by D « ksT (/? + 1) /Qitria, where (3 = 
log((a + b)/2b) (typically, (3 < 2). Thus 

5 T [\og{2 x h/b8 T )-2py 1 ' 2 

/ A \i/4 

~ X b f -J V2vr(/3+l)^ C, (18) 

Using these formulas, we have been able to fit the ex- 
perimental data on natural AFP's of two different classes: 
AFP-type I, and AFP-type III HGJ. They have rather 
different architectures: AFP I is an a-helical rod-like 
molecule which we model as a thin cylinder (b = 3.5A 
, a = 25 A), while AFP III has a more complicated 
globular structure, which we approximate as a sphere 
(b = a = 8A). Note that these lengths are derived from 
the known structures of the molecules, and therefore the 
only free parameter of our theory is x- The theoretical 
curves S T (C) are in excellent agreement with the exper- 
iments. The fact that values of the fitting parameter ( 
X oi .25 for AFP-I, and X ~ -33 for AFP-III) are phys- 
ically reasonable and close to each other gives an addi- 
tional support to our mechanism. 

Note that according to Eq. ((18(1 . the activity is 
strongly dependent on the smaller dimension of the AFP. 
This may be useful for the design of synthetic AFP's. E. 
g., by using ring-shaped molecules the S T could be in- 
creased since the effective size is set by the largest dimen- 
sion, the radius of the ring, in this case. Our results are 
consistent with the measured 6 T (C) for Antifreeze Gly- 
coproteins (AFGP). However, their molecular architec- 
ture and conformations are more complicated than those 
we have discussed, and their analysis would go beyond 
the scope of this work. 

Finally, we can compare our results with the related 
problem of the motion of a pulled elastic interface in a 
medium with static obstacles to interface motion. This 
is of interest in the description of the kinetics of domain 
walls, c harg e density waves and flux lines in supercon- 
ductors [lij • The interplay of long-range elastic coupling 
and local pinning results in a pinning-depinning transi- 
tion at a critical value of the pulling force, and near the 
transition point, the average speed of the interface goes 
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FIG. 2: (a) Simulation results for the speed of the interface as a function of time, for various values of control parameter 
fc+b 3 /uo, ranging from 0.5 x 10 -4 to 3 x 10 -4 . 5t = 0.005, \ = l/(27r) (b) Speed in the steady state vs. reduced adsorption 
rate: comparison of simulation (points) with analytical result (line) given by Eq. 1141 
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FIG. 3: Comparison of theory to experiment for maximum 
undercooling of water in the presence of AFP-Type III 
(diamonds) and two modifications of AFP-Type I (circles, 
squares) . 



continuously to zero. Our model has pinning of the inter- 
face by the AFP molecules, and coupling due to surface 
tension, but the transition is discontinuous. The differ- 
ence is that the AFP's are not stationary. Their arrival 
at the interface is controlled by the adsorption process 
and is independent of the advance of the interface. To 
emphasize this difference, we call our model 'kinetic pin- 
ning'. There should be a crossover between the static 
pinning and kinetic pinning if diffusion of the obstacles 



is taken into account. 
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